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Muon spin rotation measurements of the vortex state in vanadium: A comparative 
analysis using iterative and analytical solutions of the Ginzburg-Landau equations 
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We report muon spin rotation measurements on a single crystal of the marginal type-II supercon- 
ductor V. The measured internal magnetic field distributions are modeled assuming solutions of the 
Ginzburg-Landau (GL) equations for an ideal vortex lattice obtained using (i) an iterative procedure 
developed by E.H. Brandt, Phys. Rev. Lett. 78, 2208 (1997) and (ii) a variational method. Both 
models yield qualitatively similar results. The magnetic penetration depth A and the GL coherence 
length ^ determined from the data analysis exhibit strong field dependences, which are attributed 
to changes in the electronic structure of the vortex lattice. The zero-field extrapolated values of 
A and the GL parameter k agree well with values obtained by other experimental techniques that 
probe the Meissner state. 

PACS numbers: 74.20.De, 74.25. Ha, 74.25. Qt, 76.75.-|-i 



I. INTRODUCTION 

In order to analyze muon spin rotation (/iSR) measure- 
ments on a type-II superconductor in the vortex state, it 
is necessary to assume a theoretical model for the spa- 
tial variation of the local internal magnetic field B{r)J- 
An essential requirement of the model is that it must 
account for the finite size of the vortex cores. Thus far 
the internal magnetic field distribution n(B) measured 
by /iSR has been analyzed assuming analytical models 
for -B(r) based on London and Ginzburg-Landau (GL) 
theories. Since London theory does not account for the 
finite size of the vortex cores, a cutoff factor derived from 
GL theory must be inserted into the analytical London 
expression for B{r) to correct for the divergence of B{r) 
as r ^ 0. Unfortunately, analytical cutoff factors are 
derivable from GL theory only near the lower and upper 
critical fields, Bd and i?c2- At intermediate fields, these 
analytical cutoffs deviate substantially from the precise 
numerical calculations,^ making modified London mod- 
els inappropriate for the analysis of /iSR data. There are 
several approximate analytical expressions for B(r) that 
have been derived from the GL equationsi^i^i^ For ex- 
ample, a variational solution of the GL equations^ has 
proven to be a reliable model for analyzing /iSR measure- 
ments on VgSi (Ref.'^), NbSca (Rei.^) and YBaaCugOy-a 
(Refi^*i£). However, this analytical GL model is strictly 
valid only at low reduced fields b = B/Bc2 and large val- 
ues of the GL parameter k. Thus, often used analyti- 
cal models for B{r) have limited validity and can devi- 
ate substantially from the numerical solutions of the GL 
equations. 

Brandt has developed an iterative method for solving 
the GL equations that accurately determines B{r) for 
arbitrary b, k, and vortex-lattice symmetry. ^^ Thus far 
this iteration method has not been applied to the analy- 
sis of /iSR measurements of n(B) in the vortex state. As 
a first test of this method we have chosen to study the 
marginal type-II superconductor vanadium (V). This rig- 



orous analysis method is expected to be required for V, 
whose low value of k falls outside the range of validity of 
the analytical models. In addition, the low value of Bc2 
(« 0.45 T) gives us experimental access to reduced fields 
which are beyond the range of validity of the analytical 
model. 

The paper is organized as follows: Theoretical models 
for B{r) are described in Sec. II. The experimental pro- 
cedures are described in Sec. III. Measurements in zero 
external magnetic field are presented in Sec. IV. Mea- 
surements in the vortex state are described in Sec. V 
and concluding remarks are given in Sec. VI. 



II. THEORETICAL MODELS 

A. Iterative GL solution 

Here we briefly outline the iteration method presented 
in Ref.^^, and correct some typographical errors con- 
tained therein. The GL equations are written in terms 
of the real order parameter lo, the local magnetic field 
B, and the supervelocity Q, which are expressed as the 
Fourier series 



K 



COS K • r) 



B{r) = B + ^bKCOsK-r, 
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Q(r) = Q^(r) + ^ ^k^^^ sinK • r 
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(1) 
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where ok and &k are Fourier coefficients, w(r) — |?/>(r)p, 
B is the average internal field, and '(/'(r) is the complex 
GL order parameter. The "tail" of the position vector 



r = {x,y) is at the vortex center (0,0). The local mag- 
netic field B is given in units of V^Bc (where Be is the 
thermodynamic critical field) and all length scales are in 
units of the magnetic penetration depth A. Q^(r) is the 
supcrvclocity obtained from Abrikosov's solution of the 
GL equations near Bc2 
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where ujAiic) is calculated from Eq. (Q using 



a^ = -{-!)" cxp(-7ri^%/3). (5) 

^ assuming a hexagonal vortex 



Here i/ = m + mn + n 

lattice with vortex positions given by 
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(mxi +nx2,ny2), 



(6) 



where m and n are integers, xi is the intervortex spac- 
ing, X2 = a;i/2 and j/2 = a;i\/3/2. The spatial field pro- 
file B{r) was calculated at approximately 950 locations 
within one-quarter of the rectangular unit cell shown in 
Fig. n^a) . The reciprocal lattice vectors used in the cal- 
cination of B{r) are given by 

2tt 
K = K„in = -^{my2,nxi -I- 771x2), (7) 

where S is the unit cell area. The Kmn vectors are re- 
stricted to those indicated in Fig.^b), corresponding to 
— 16 > 777 < 16 and — 16 > ri < 16 within a semicircle with 
Ky > (but excluding vectors with Kx>0 and Ky = 0). 
It was found that the calculated field distribution did 
not change significantly if the summation was extended 
to values of |77| and |77i| greater than 16. 

The Fourier coefficients ok and 6k are calculated from 
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FIG. 1: (a) The vortex positions in the hexagonal lattice (O) 
and the rectangular unit cell (dotted lines). Due to sym- 
metry, the theoretical probability field distribution n{B) was 
calculated at approximately 950 locations (•) within just one- 
quarter of the rectangular unit cell. For the sake of clarity, 
the sampled positions are shown at 1/4 of the density used in 
the calculations, (b) The points in the reciprocal space lattice 
(o) that are included in the K sums of Eqs. Q-Q, Eq. IIH 
and Eq. ifT^ . 
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6k = 



-2 < [ujB + u){B - B) + p] cos K • r > 



K^ 



is the spatial average of lo 
- QxjT- We note that the 



where g — {Vtuy /{in'^uj) 

and p= (Vw X Q)z = Qy-Q^ - ^x^^ 

definitions of p and the Cj{B — B) term in Eq. Hl()|l are 

incorrectly written in the original article, ^2 i^^it have been 

corrected here. As explained in Refiii, solutions to the 

GL equations are acquired by first iterating Eqs. 1^, © 

and a few times to relax w and then iterating Eqs. 

(Uni, (HJ, Q, ®, ©, ® and again ^^ etc. . . to relax 

B. 



B. Comparison with other models 

Here we compare the results of the above iteration 
method for B{r) to the widely used modified London 
and analytical GL models. The local magnetic field at 



(10) position r={x,y) in the modified London modcl^^ is 



S(r) = B^. 



K 



iK-r^-K^f/2 

TTA2^2 ■ 



(11) 



Although this model is considered applicable for reduced 
fields b = B/Bc2 < 0.25 and k > 2, the Gaussian cutoff 
factor exp(— iir2^2y'2) introduced to account for the log- 
arithmic divergence of B{r) at the center of the vortex is 
not strictly valid. 2 

The approximate analytical solution of the GL equa- 
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FIG. 2: (Color online) The spatial field profile B{r) along 
the straight line connecting nearest-neighbor vortices, for an 
hexagonal vortex lattice, k = 25, B = 3 kG and A =1000 A. 
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FIG. 3: (Color online) The spatial field profile B(r) along 
the straight line connecting nearest-neighbor vortices, for an 
hexagonal vortex lattice, k.— \/\/2, 5 = 753 G and A =150 A. 



tions for B{y) \Sm 

B{y) = 5(1 - b 



4,^e-''^-'^uKi{u) 



K 
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(12) 



where 



^2^^K^{l + b'^)[l-2b{l-bf 



(13) 



Ki is a modified Bessel function and ^ is the GL coher- 
ence length. This analytical GL model is a reasonable 
approximation for 6^1 and k3>1. 

Figures [21 El and ^ show comparisons between the so- 
lutions for B{r) from the three different models, plot- 
ted along the straight line connecting nearest-neighbor 
vortices. The parameters used to generate the curves 
in Fig. 121 are characteristic of the high-K superconduc- 
tor VsSi (Reffi). While there is good agreement between 
the iterative and analytical solutions of the GL equations, 
the modified London model deviates substantially in the 
region of the vortex cores. Figure shows that the mod- 
ified London and analytical GL models completely break 
down for k = 1/\/2, the limit of type-II superconductiv- 
ity. For example, in the case of the analytical GL model, 
there is actually a region between the vortices where the 
local field changes direction. In Fig. 01 plots of B(r) are 
shown for a set of parameters obtained from fiSK mea- 
surements on the low-K superconductor V (see Sec. V). 
The value k = 5.3 is rather large for V, but as we explain 
in Sec. V, k is really an "effective" fit parameter influ- 
enced by the electronic structure of the vortex lattice. 

In Fig. 01 the analytical GL solution deviates from the 
iterative GL solution both near and far from the vortex 
core regions, whereas the modified London model fails 
only in the region of the vortex cores. 
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FIG. 4: (Color online) The spatial field profile B{r) along 
the straight line connecting nearest-neighbor vortices, for an 
hexagonal vortex lattice, k = 5.3, _B = 1.6 kG and A =1042 A. 



III. EXPERIMENTAL DETAILS 

The single crystal of the 1ow-k superconductor V mea- 
sured in the present study was purchased from Goodfel- 
low Cambridge Ltd.'^^^ The crystal is a disk, 13 mm in di- 
ameter by 0.35 mm thick, with the (111) crystallographic 
direction perpendicular to the plane of the disk. Magne- 
tization measurements indicate that the crystal has a su- 
perconducting transition temperature of Tc = 5.2 K, and 
an upper critical field of Hc2 ~ 4.2 kOe. This value of Hc2 
corresponds to a BCS coherence length of ^o ~ 280 A. 
A four-probe potcntiometric ac resistivity measurement 
yielded p — 0.7 /iO-cni just above Tc, which is 31 times 
smaller than p at rrom temperature. Using the carrier 
concentration n « 9 x 10^* m~'^ obtained from Hall resis- 
tance measurements^ a free electron theorjiiS. calcula- 
tion of the mean free path yields I = hkp/ pne^ w 900 A, 
where kp = (37r^n)^''^ is the Fermi wavenumber and e 



is the electronic charge. Thus our sample is in the clean 
limit with //^q ~ 3. Neutron scattering measurements 
performed on our V single crystal showed no evidence of 
an intermediate mixed state {i.e. a mixture of Meissner 
and vortex-lattice phases) lii 

The /iSR measurements were carried out on the M15 
beamline at the Tri-University Meson Facility (TRI- 
UMF), Vancouver, Canada, using a dilution refrigera- 
tor to cool the sample. Measurements of the vortex 
state were done under field-cooled conditions in a "trans- 
verse field" (TF) geometry, in which the magnetic field 
was apphed along the z-axis parallel to the (111) di- 
rection of the crystal, but perpendicular to the initial 
muon spin polarization Px{0) (which defines the x-axis). 
Each measurement was done by implanting approxi- 
mately 2x10^' spin-polarized muons one at a time into the 
crystal, where their spins precess around the local mag- 
netic field B{r) at the Larmor frequency lo = ^f^B, where 
7^ = 0.0852 /^s^^ G~^ is the muon gyromagnetic ratio. 
The muons stop randomly on the length scale of the vor- 
tex lattice, and hence evenly sample B(y). The /^SR sig- 
nal obtained by the detection of the decay positrons from 
an ensemble of muons implanted into the single crystal 
is given by 



A{t) = aaPx{t) ^ 



(14) 



where A{t) is called the /iSR "asymmetry" spectrum, uq 
is the asymmetry maximum, and Px{t) is the time evo- 
lution of the muon spin polarization 



/>OC 

Px (t) = / n{B) cos(7^Si + <j))dB . 
Jo 

Here </> is a phase constant, and 

n{B') = {6[B' - B{r)]) , 



(15) 



(16) 



is the probability of finding a local field B in the z- 
direction at an arbitrary position r in the x-y plane. Fur- 
ther details on this application of the /iSR technique are 
found in Refi. 
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FIG. 5: (a) Representative asymmetry spectra (symbols) ac- 
quired in V in zero external magnetic field, and fits (solid 
lines) to the numerical dynamic Gaussian Kubo-Toyabe func- 
tion, (b) The temperature dependence of the extracted muon 
hop rate v (circles). The solid line is a fit to ^{T > 15 K) 
using an Arrhenius function. 



IV. ZERO-FIELD MEASUREMENTS 

Figure[SJa) shows asymmetry spectra for our V sample 
in zero external magnetic field. These spectra contain 
a 7 % time-independent background contribution from 
muons stopping outside the sample. The signal com- 
ing from muons stopping inside the sample is well de- 
scribed by a numerical dynamic Gaussian Kubo-Toyabe 
function.^* This function is characterized by a relaxation 
rate A corresponding to the width of the internal mag- 
netic field distribution experienced by the muons, and a 
parameter v corresponding to the hop rate of the muons 
in the sample. As shown in Fig. |3Jb), the muon hop 
rate in our V crystal decreases with decreasing temper- 
ature to i/wO.2 ^s-i at T = 15 K (l/T w 0.07 K^^). 



At r > 15 K, the data are well described by the classi- 
cal Arrhenius law for thermally activated motion in the 
presence of potential barriers^- 



i/Qexpi-Ea/kBT), 



(17) 



where ks is Boltzmann's constant, Ea is the activation 
energy for thermally assisted muon hopping, and vq is a 
constant. Fitting the T > 15 K data using this expression 
yields vq = 1.72 ^s^^ and Ea = 4.5 meV. Below T = 15 K 
there is perhaps a slow increase in the hop rate v, which 
we speculate is due to quantum mechanical tunneling as 
observed in other metalsiiS 

Assuming that the muon occupies an interstitial site 
of tetrahedral symmetry in the V crystal lattice, we can 



calculate the muon diffusivity Dn from the expressioi 



,18 



D, = -24, 



(18) 



where 

D^ « 9.7 X 10 



3.02 A is the lattice constant. This gives 
1^ mVs at T = 2.5 K. Brandt and 
Seeger performed a thorough theoretical study of the ef- 
fect of muon diffusion on //SR lincshapes in the vortex 
stateiSS They found that muon diffusion causes signifi- 
cant smearing of the sharp features of n{B) for values 
of D^ greater than « 10~^7^|M|d^, where M is the 
sample magnetization and d is the intervortex spacing. 
Our measured muon diffusivity is several orders of mag- 
nitude smaller than this. At if = 1.6 kOe, for example, 
we have 7^|M|d^ w 2 x 10~® m^/s, which means that 
^^(w 9.7 X IQ-^^ mVs) = 5 X lQ-^-f^\M\(P. Thus muon 
diffusion has a negligible effect on the /iSR linehapes mea- 
sured here. 
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FIG. 6: (Color online) Time evolution of the muon spin polar- 
ization (cirlces) in V at H = 1.Q kOe and r = 2.6 K, and a fit 
(solid curve) assuming the solution of B{r) from the iterative 
GL method. 



V. MEASUREMENTS IN THE VORTEX STATE 
A. Comparison of fits 

To fit the /iSR signals in the vortex state, the field dis- 
tribution n{B) contained in the theoretical polarization 
function P^ [t] of Eq. H15|l was generated from one of the 
three theoretical models for B{t) described in Sec. II. For 
all cases, an ideal hexagonal vortex lattice was assumed. 
In addition, P^ (t) was multiplied by a Gaussian function 
exp(— cr^t^/2), which is equivalent to convoluting n{B) 
with the Gaussian exp[— (7^5/(7)^/2]. This accounts for 
disorder in the vortex lattice^ and the static local-field 
inhomogeneity created by the large ^^V nuclear dipole 
moments. An additional Gaussian depolarization func- 
tion was added to Eq. l(n|) to account for approximately 
20 % of the signal arising from muons that stopped out- 
side the sample. 

A typical asymmetry spectrum at iJ = 1.6 kOe and 
T = 2.6 K is displayed in Fig.O The solid curve through 
the data is a fit assuming the solution of B{r) obtained 
from the iterative GL method. The parameter values 
obtained from this fit were used to calculate the spatial 
field profiles shown in Fig. 01 Fourier transforms of typi- 
cal time domain signals and fits to both the iterative and 
analytical models are shown in Figs. [T] and |H1 From the 
Fourier transforms one can see that both fits capture the 
main features of the /iSR lineshape and are of high qual- 
ity. In particular, for the data in Fig.[7|the ratio of x^ to 
the number of degrees of freedom (NDF) is comparable, 
being 1.26 for the fit to the iterative GL solution and 
1.30 for the fit to the analytical GL model. The values 
of A, ^ and a extracted from the two models differed by 
9 %, 8 % and 13 %, respectively. We note that despite 
returning significantly different parameter values in some 
cases, fits with both models resulted in similar values of 
X^/NDF for all of the data presented in this article. 



Interestingly, the results from the two models are in 
slightly better agreement at the lowest temperatures and 
magnetic fields. For example, at H =1.2 kOe and T = 
0.02 K, the differences in A, ^ and a are 7 %, 8 % and 
2 %, respectively. On the other hand, at if = 2.9 kOe 
and T = 0.02 K, A, ^ and a differ by 10 %, 9 % and 3 %, 
respectively. Even so, the quality of the fits is about the 
same for both models. This is evident from the Fourier 
transforms shown in Fig. |S1 At high reduced field b, the 
values of A and ^ obtained from the iterative GL method 
are likely to be more accurate, since at these high fields 
the analytical GL model is being applied outside its range 
of validity. 

In the next section, complete results for fits using both 
the iterative and analytical GL models are presented. 
There we show that despite differences in the absolute 
values of A and £,, fits to both models yield similar tem- 
perature and magnetic field dependences for these length 
scales. In particular, we show that the value of A extrap- 
olated to zero field is, within experimental uncertainty, 
identical for both models, and furthermore agrees well 
with Meissner state measurements on V using other ex- 
perimental techniques. 



B. Temperature dependences of A and ^ 

Fourier transforms of the muon spin precession signal 
from V at ii = 1.6 kOe and temperatures below Tc are 
shown in Fig. [^ Magnetization measurements indicate 
that Tc = 3.65 K at ii = 1.6 kOe. As the temperature 
is lowered, the /.tSR line shape broadens and the ampli- 
tude of the high-field "tail" decreases. While the high- 
field cutoff is less obvious at T = 0.02 K, we note that 
the "true" cutoff in n{B) is smeared out by the Fourier 
transform^ In fact the fits in the time domain are quite 
sensitive to the high-field tail, yielding finite values for ^, 
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FIG. 7: (Color online) Fourier transforms of the muon spin 
precession signal in V at H = 1.6 kOe and T — 2.6 K (circles), 
the fit using the iterative GL solution (blue curve) and the fit 
to the analytical GL model (red curve). The peak at 1600 G 
is the background signal originating from muons that missed 
the sample. 




2700 2800 



2900 3000 
B(G) 



3100 3200 



FIG. 8: (Color online) Fourier transforms of the muon spin 
precession signal in V at H = 2.9 kOe and T = 0.02 K (circles), 
the fit using the iterative GL solution (blue curve) and the 
fit to the analytical GL model (red curve). The peak near 
2900 G is the background signal originating from muons that 
missed the sample. 



even at T = 0.02 K. 

Figures HTJl and ITTI show the temperature dependences 
of 1/A^, a and ^ determined from our fits of the muon 
spin precession signals at H ^ 1.6 kOe, assuming solu- 
tions for B{r) from both the iterative and analytical GL 
methods. Despite the differences in absolute values of 
1/A^, both data sets in Fig. ^1 are well described by 
BCS weak-coupling 1/\^{T) curves^^ for Tc = 3.65 K. 
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FIG. 9: Fourier transforms of the muon spin precession signal 
in V at H = 1.6 kOe and T<Tc{H) = 3.65 K. 



The inset in Fig. ^| shows that both models yield sim- 
ilar values for the additional broadening parameter a. 
As A becomes longer with increasing temperature, there 
is a greater overlap of the vortices and a corresponding 
reduction in the pinning-induced disorder of the vortex 
lattice. This is because the energy associated with the 
interaction between vortex lines depends on A.«3 For this 
reason a{T) roughly follows 1/A^ in Fig. IIUI 

The temperature dependence of the GL coherence 
length ^ is shown in Fig. ^2 ■C is a measure of the 
vortex core size. Recently, we have demonstrated from 
fiSH and thermal conductivity measurements on BCS 
superconductors^iS that the core size is dependent on 
the degree of localization of the quasiparticle bound core 
states. Thermal depopulation of the more spatially ex- 
tended high-energy core states results in a shrinking 
of the core size with decreasing temperature. This is 
the so-called "Kramer-Pesch effect",^'' which has pre- 
viously been observed in NbSe2 by /^SR^'^^ and shown 
to be dependent on magnetic fieldi2& In a clean BCS 
type-II superconductor the core size of an isolated vor- 
tex is expected to be temperature independent below 
T^ Tc/kp^Q, where kp is the Fermi wave number and 
^0 is the BCS coherence length. We see in Fig. E] 
that ^(T) obtained from the fits to both models dis- 
plays the Kramer-Pesch effect, with ^ saturating be- 
low T « 1 K. Using the free-electron expression kp = 
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FIG. 10: Temperature dependence of 1/A^ in V at // = 
1.6 kOe, determined from fits using the iterative and analyt- 
ical GL models. The solid and dotted curves are theoretical 
BGS weak-coupling predictions for rc = 3.65 K. 
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FIG. 11: Temperature dependence of ^ in V at H = 1.6 kOe, 
determined from fits using the iterative and anlytical GL mod- 
els. 



(37r2n)i/3 (Ref.i6) and n w 9 x lO^s m^^ from Hall resis- 
tance measurements jiS. we obtain kp ~ 1.4 x 10^° m~^. 
Assuming the value of the superconducting coherence 
length ^0 ~ 280 A estimated from the extrapolated zero- 
temperature value of Hc2, the core size in our V crys- 
tal is therefore theoretically expected to saturate below 
T « 10^^ K. The premature saturation of the core size 
observed in Fig. II II could result from quasiparticle scat- 
tering by nonmagnetic impurities,'^- but this is unlikely 
given that our sample is in the clean limit. It is important 
to note that theoretical predictions only exist for isolated 
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FIG. 12: Fourier transforms of the muon spin precession sig 
nal in V at T = 0.02 K and H<Hc2- 



vorticesi2^*2LS£ In a lattice, the core states of nearest- 
neighbor vortices overlap to some degree, and this is 
likely the reason why the strength of the Kramer-Pesch 
effect observed by /iSR weakens with increasing fieldi2& 
The delocalization of core states due to vortex- vortex in- 
teractions, and the corresponding reduction in the core 
size, also explains why the low-temperature value of ^ in 
Fig. E] is much smaller than <^o- 



C. Magnetic field dependences of A and ^ 

In Fig. 1121 Fourier transforms of the muon spin pre- 
cession signal in V at T = 0.02 K are shown for different 
applied magnetic fields H < if c2 • The changes in the fiSR 
line shape as a function of H are similar to that previ- 
ously observed in NbSe2 (Ref j^^) , and result directly from 
the change in vortex density. Increasing the vortex den- 
sity reduces the internal magnetic field inhomogeneity 
and increases the degree of overlap of the wave functions 
of the core states of neighboring vortices. 

The magnetic field dependence of A in V, as deter- 
mined by both the analytical and iterative GL models, 
is plotted in Fig. 1131 along with our previously published 
data for VgSi (Ref.^) and NbSca (Ref.«). A from both 
models is seen to increase with field, although use of the 
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FIG. 13: Magnetic field dependence of A determined from 
fiSR measurements on V, VaSi (Ref.'') and NbSe2 (Refp&). 
The straight lines are linear extrapolations of the data to zero 
field. Inset: Blow up of the low-field data for VsSi and NbSe2. 



iterative model results in a slightly stronger field depen- 
dence. The value of A determined by /iSR depends on 
the radial decay of B{r) outside the vortex cores. Since 
the spatial field profile around a vortex core can be sig- 
nificantly modifed by the delocalization of bound core 
states, the measured value of A may be strongly depen- 
dent on field. We stress that when this is the case, X{H) 
is an "effective" length scale, which in the fits absorbs 
changes in B{r) due to changes in the electronic struc- 
ture of the vortex lattice. This dominates over the weak 
field dependence of A expected for an isolated vortex in 
an s-wave superconductor.-*? We note that the difference 
in slope of A vs iJ between the two models in Fig. E| 
suggests that the exact details of how these changes in 
electronic structure are absorbed by X{H) is model de- 
pendent. To compare with measurements of A by other 
experimental techniques, we have extrapolated the data 
for X{H) to H^>0 kOe. The zero-field extrapolated value 
of A in V is 376. 3± 22.4 A using the iterative GL model 
and 375. 9± 17.0 A using the analytical GL model. Mag- 
netization measurements'^^-'^^'''^ have determined that A 
is in the range 374 A to 398 A, in excellent agreement 
with our results. We see that either model for B{r) can 
be used to extract a reliable measure of the zero-field 
magnetic penetration depth. 
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In VsSi, where the bound core states are highly lo- 
calized at low fields^ A is weakly dependent on field 



FIG. 14: Magnetic field dependence of ^ in V at T = 0.02 K 
as determined from fits using the iterative and analytical GL 
models. 



below H K,Q.2Hc2 (RefZ). The zero-field extrapolation 
shown in the inset of Fig. [H yields A = 1080 ± 17 A, 
in good agreement with the value A= 1050 A determined 
from Hc2{T) measurements in Ref.'''^. Likewise, the zero- 
field extrapolated value A = 1249±31 A for NbSe2 agrees 
well with the value of 1240 A obtained in Refi^ from 
magnetization measurements. Thus /iSR can be used 
for accurate measurements of the absolute value of the 
magnetic penetration depth in type-II superconductors, 
provided there is sufficient data to permit an accurate 
extrapolation to zero field. This works even in the case 
of an unconventional superconductor. Recently, it was 
shown that zero-field extrapolated values of A obtained 
from /^SR measurements on the high temperature super- 
conductor YBa2Cu306+a: agree well with values obtained 
from accurate electron spin resonance measurements in 
the Meissner phasei^ We note that the linear extrapo- 
lations of the data in Fig. E| are continuous through the 
Meissner phase — which occurs in V below iJci ~0.25iJc2- 
The magnetic field dependence of ^ at T = 0.02 K is 
plotted in Fig. ^1 Although the analytical model yields 
larger values of ^, both models display a similar field de- 
pendence over the entire field range. Immediately above 
Hci, the vortex core size (~ ^) shrinks with increasing 
field and saturates above H « 2.4 kOe. (We note that 
analysis of a recent Andreev reflection spectroscopy study 
of niobium revealed a similar trend for £,{H) over the 
same range of reduced fields 0.3 < b < 0.7.'^^) We at- 
tribute this behavior to an increase in the overlap of the 
core states of nearest-neighbor vortices;^ as was found 
to be the case in VaSi and NbSe2P^ Although Kogan 
and ZhelezinaiS, have also successfully modeled the field 
dependence of the core size in VsSi and NbSe2 by weak- 
coupling BCS theory, their calculations assume a large 
GL parameter k — A/^ and hence are not applicable 
to V. In an isotropic s-wave superconductor, the delo- 
calization of core states is predicted to be significant 
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FIG. 15: Magnetic field dependence of k = A/^ in V at T = 
0.02 K as determined from the iterative and analytical GL 
models. The solid (dotted) line is a linear fit to the data 
obtained from the iterative (analytical) GL model (see main 
text). The fit to the data obtained from the analytical model 
excludes the points at H = 2.65 kOe and 2.9 kOe. 



at fields above B* w Bc2/i, although the value of this 
crossover field is reduced somewhat by anisotropy.**^ Spe- 
cific hea1^- and ultrasonic attenuation-^?- measurements 
suggest that the anisotropy of the superconducting en- 
ergy gap in V is approximately 10 %. According to the 
calculations of Ref.*^^, significant delocalization of the 
core states and a reduction in the core size should oc- 
cur above B* « 1.3 kG. Given the uncertainty in the 
values of Bc2{T -^ 0) and the anisotropy, the observed 
shrinking of the vortex cores at fields H > 1.2 kOe seems 
reasonable. 

Finally, we present the magnetic field dependence of 
the GL parameter k = X/^ in Fig. E| The increase 
in K with field is due to the field dependences of A and 
£_. The straight lines in Fig 1151 show that k is roughly 
a linear function of H. Significant deviations from this 
behavior occur for the values of k at the two highest 



fields determined from fits to the analytical GL model. 
This is perhaps due to the fact that the analytical model 
breaks down at high reduced fields b. It was found that 
both scatter and uncertainty in the values of A and ^ were 
considerably reduced by fixing n at each field to lie on the 
straight line fits shown in Fig. E| The ^{H) and X{H) 
data in Figs. 1131 and 1141 were obtained in this way. The 
zero-field extrapolated values of k determined from the 
fits shown in Fig.^jare 1.26±1.01 for the iterative model 
and 1.30 ± 0.64 for the analytical model. We note that 
both of these values agree well with the value k = 1.34 
calculated from our estimated values of X{H = 0) and ^o- 
Also, within experimental uncertainty, both extrapolated 
values of k are comparable with that obtained by other 
experimental methods for samples of similar purity^ 
VI. SUMMARY AND CONCLUSIONS 

We have analyzed /LtSR measurements of the internal 
magnetic field distribution in the vortex state of the low- 
K type-II superconductor V using both Brandt's iterative 
GL methodii and the analytical GL model of Ref.''. Sur- 
prisingly, the two models produce qualitatively similar 
results for both the temperature and field dependences 
of A and ^. In particular, fits to each model yield low- 
temperature, zero-field extrapolated values of A and k 
that agree with previous measurements of these quanti- 
ties by other techniques. We find that the largest dif- 
ference between the results using these models occurs at 
high fields, where the analytical GL model is being ap- 
plied outside its range of validity. The observed field 
dependences of A and ^ in V are likely due to the de- 
localization of quasiparticle core states, as has already 
been established in other conventional superconductors. 
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